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H ! Abstract 

We consider the small mass asym.ptotics (Sm.oluchowski-Kramers approxima- 
tion) for the Langevin equation with a variable friction coefficient. The limit of 
the solution in the classical sense does not exist in this case. We study a mod- 
ification of the Smoluchowski-Kramers approximation. Some applications of the 
n ', ' Smoluchowski-Kramers approximation to problems with fast oscillating or discon- 

f-H ' tinuous coefhcients are considered. 
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J5 ■ 1 Introduction 

p 

«^ ' The Langevin equation 

o 

(N ; /.qI' = 6(9^ - A^r + c7{q>;)Wt , q^ = q G M" , q^ = p G M" , (LI) 

describes the motion of a particle of mass fi in a force field b(q), q G M", subjected 
^ \ to random fluctuations and to a friction proportional to the velocity. Here Wt is the 

C^ , standard Wiener process in M", A > is the friction coefficient. The vector held h{q) 

and the matrix function o"(q) are assumed to be continuously differentiable and bounded 
together with their first derivatives. The matrix a{q) = {aij{q)) = cr{q)a*{q) is assumed 
to be non-degenerate. 

Put p^ = (jf. Then (1.1) can be written as a first order system: 
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The diffusion process (pf , q^) = X^ in M is governed by the generator L: 
Note that, since functions q^ are continuously differentiable with probabihty one, 

JO 

This allows to consider equations (1.2) for each trajectory VFj individually, and 
there is no necessity in the introduction of a stochastic integral. In particular, if (1.2) 
is considered as a stochastic differential equation, stochastic integrals in the Ito and in 

the Stratonovich sense coincide: / a^q^^^dW g = / cr{q^) o dW g- 

Jo Jo 

It is assumed usually that the friction coefficient A is constant. Under this assump- 
tion, one can prove that q^ converges in probability as ^u J, uniformly on each finite 
time interval [0,T] to an n-dimensional diffusion process qf for any k,T > and any 
p{; = p G M" fixed. 



limP max Iq'r — qAmd > k ] = . 

f^iO \0<t<T'^ ''"^ J 

Here q^ is the solution of equation 

Qt = \biQt) + l^iQtWt , Qo = g^ = 9 e IR" • (1-3) 

The stochastic term in (1.3) should be understood in the Ito sense. 

The approximation of q^ by q^ for < ^ << 1 is called the Smoluchowski-Kramers 
approximation. This is the main justification for replacement of the second order equa- 
tion (1.1) by the first order equation (1.3). The price for such a simplification, in 
particular, consists of certain non-universality of equation (1.3): The white noise in 
(1.1) is an idealization of a more regular stochastic process Wf- with correlation radius 
5 « 1 converging to Wt as 5 ]. 0. Let q^' be the solution of equation (1.1) with Wt 
replaced by W^ . Then limit of q^' as /i, (5 J, depends on the relation between fi and 
6. Say, if first (5 | and then fi iO, the stochastic integral in (1.3) should be understood 
in the Ito sense; if first /x | and then 5 | 0, q^' converges to the solution of (1.3) with 
stochastic integral in the Stratonovich sense. (See, for instance, [5].) 

Consider now the case of a variable friction coefficient A = A(q). We assume that 
A(g) has continuous bounded derivatives and < Aq < A(q) < A < oo. It turns out, 
as we will see in the next section, that in this case the solution q^ of (1.1) does not 
converge, in general, to the solution of (1.3) with A = A(q), so that the Smoluchowski- 
Kramers approximation should be modified. In order to do this, we consider equation 
(1.1) with Wt replaced by W^ described above: 



m'/ = h{q^/) - X{q^'')q^'' + ^{q^'^Wt , qf = q , qi;'' = P ■ (1.4) 

We prove that after such a regularization, the solution of (1.4) has a hmit q^ as 
/i 4, 0, and Qj is the unique solution of the equation obtained from (1.4) as ;U = 0: 

Qt = T^Mqf) + -±^a{qf)wt , q'o = q ■ (1.5) 

Now we can take 5 4, in (1.5). As the result we get the equation 

Qt = YJ^rHQt) + w^Tf^(Qt) °Wt , qo = q , (1.6) 

where the stochastic term should be understood in the Stratonovich sense. So the 
regularization leads to a modified Smoluchowski-Kramers equation (1.6). We prove this 
in Section 3. 

Some applications of the Smoluchowski-Kramers approximation are considered in 
Sections 4 and 5: the case of fast oscillating in the space variable, periodic or stochastic, 
friction coefficient is studied; gluing condition at the discontinuity points of the friction 
coefficient are considered. In the last Section 6, we briefly consider some remarks and 
generalizations. 

Notations. We use | • \^d to denote the standard Euclidean norm in M"^. When 
d = 1 we set I • 1^1 = I • |. For a vector-valued function f{x) = {fi{x), ..., fd{x)), x G R^, 
we set 11/ II 00 = max ||/j||oo = max sup |/j(a;)|. All the vectors are marked with either 

l<i<d 'i-<i<d^^^d 

bold letters or with an arrow on it. 



2 Some estimates. The classical Smoluchowski-Kramers 
approximation does not work for variable friction coef- 
ficients 

We consider the following system 

^qf = b(qn - A(qr)gr + Wt, q^o=Q^^'^Qo=P^^'- (2.1) 

Here 00 > A > A(») > Aq > is a function of q^. We assume that function A(«) and 
the vector field 6(») are continuously differentiable and bounded together with their 
first derivatives. The process Wt is the standard Wiener process in M . For simplicity 
of calculations we consider here the case when the diffusion matrix a(») is the identity 



(compare with (1.1)). The case of general diffusion matrix can be considered in a similar 
way and we will briefly mention it in Section 6. 

Let p^ = q^, we have, that (2.1) is equivalent to the system 

Pt = -M€) — j^Pt + -Wt . 

Then 



^ ^eUo^^i^)''-'p'^) ^ gi/o*A(9!^)rfs Tm + }-X(qi^)pA = e^/o ^('?^)'^« f-b(gn + -ii^t 



dt \ J \ fi J \fi fi 

and 

eTi /o A(<?^)rfSpM _ p = i /"* e^ /o ^('?^)'^'-b(q^)ds + i /"* cm /o ^(''r)'^^ dW s . (2.3) 

Ai Jo Ai Jo 

For notational convenience we introduce the function A[fi,t) = / X{q'^)ds. It is 

Jo 
clear that tA > A(^, t) > tAo. Using (2.3) we have, that 



^ Jo M Jo 

Therefore we have 



qt=<i+ [ p^ds 

J,o 



= q+p l' e-Ti^^^^'^ds +- f e-^^(^'^) ( f e^''^^'^''h{q^,)dr\ ds+ ^ ^ 

Jo f^ Jo \Jo J (2.4) 

= 9 + «(/") + /9(/") + 7(/") • 
Here q;(/u),/3(^),7(/u) are three (vector) functions in the right hand side of (2.4): 



-j:A{^^,s) 



ds 



a{iJ,) = p e m' 
Jo 

/3(^) = 1 re-7r^('^'^) ( r e^^^^''''h{q^)dA ds , 



In the following we will use the relation 

d_ 

dt \ / /i 

We will also use the estimates 



V / /i dt /i ^ * ^ ' 



(2.5) 



— -(1-e M 
cA 



=-S^(M,^), 



A* 



Ai 



e i^ ds < I e i^ 'ds < / e i^ ds = — — (1 — e '')<,, 
Jo Jo cXq cXq 

(2.6) 



cA 



eA(t-s) 



e M ds < I e t^ 







-^(A(t^,t)-A(p,s)) 



ds < 







< e M ds = -V-(l-e '^ ) < -^ . 
Jo cAo cAo 



(2.7) 



Here c is a positive constant. 



We get in this section some bounds for a(/x), /3(/i), 7(/u) which show, in particular, 
that the classical Smoluchowski-Kramers approximation does not hold in the case of 
variable friction. These bounds also will be used to obtain a modified Smoluchowski- 
Kramers approximation. 

2.1. Estimates of a{fi). 



We have, by (2.5), 

"(^) =P e' 



.-^^(M,s) 



ds 



-PH 



-A{^l,t) 



1 



KqD A(g) 7o ' ' ^A(q^) 



Let Rctip) = A* 

|-Ra(/u)| — )■ as /x I 0. 
Let 



We have 



e M 



A(M,t) 



1 



A(gr) A(g) 



^ 



. It is easy to estimate | Ra (^) | < t— • Therefore 

Ao 



{I) 



Kqs) 



{I) 



e---^^''''^^^^VX{q^)-p^ds 



e '^ 



p + 1 rei^('^'^)/,r«M 






(/i) + ih) + (/s) . 



Here 



ih 



-p • e M 

/o 



* lA(M,s)VA(q^ 



A2(q: 



ds 



(I,) = -i / e-t^(^'^)-^VA(g^) • ( r e^^^^^Mq'r)dr 1 ds , 



(/3 



1 '■* 



A2(q^^ 



-— e 1^ 



iAi^,s)_ 



XHq: 



-VX{q: 



e-.^^'^^'UWr 1 ds . 



We can derive, using (2.6) and (2.7), that 



Since 



(^i)l < — 79 — IpI: 



\2 
^0 



e-^^^'ds < 



\ih)\< 



IVAI 



'0 

1 '•* 



llVAIIc 



JJ_ 
'2An 



^ l|VA||oO||.|| 

-^0 



--A{n,s) 



ei^(^''-). 







IVA 



^ Jo 



dr ds 



< 



< 



VA 



T2 — l|o||oo- / -^(1-e '^ )e ^^ ds 
Aq M jo -^0 



VA 



Ag 



loo / e M ds 
'0 

l^llooT- • 
'^O 



1(^3)1 < 



|VA||ool 
Aq /^ 



/■ e-^-^(^'^) (" re-^^(^'^)e-^^('^'^)+^^('^'^)d\^.] ds 



we could estimate, by Cauchy-Schwarz inequality and (2.6), (2.7), that 



E|(/3)|'< 



IVAI 



^E 



< 



< 



< 



< 



< 



X2 / ,,2 

VA|U\^ 1 ' '■' 



\~h^M ( /'%-^A(^,s)^^iA(^,s)+iA(^,r)^ , ^^ 



X2 

VAII 



/^ Vio 



E / e—y.^y^^'>ds 







--A(m,s) 

e '^ 



^-lA(^,.)+lA(M,r)^^^ 



ds 



X2 

VAII 



/^ VJo 



e M (is 



X2 
VA 



/^ VJo 



— 77 I / e M (is 



OO \ L I I U_ 



e A* E 



e M 



^-lA(M + lA(^,r)^^ 



(is 



Ee" 



■|A(M,s)+|A(M,r)^ 



dr I ds 



e ^ (is 










e ^ 



-+ 



e ^ M (ir ds 



Aq ) /^^ Ao 

VA||oo \ /i 

Ag ) 2A^ 







2A, 







Combining these estimates we see that E|(/)p — )■ as ^ J, 0. 
So E|q!(^)|^^ — )• as /U 4, for any |p|]Rd < 00. 

2.2. Estimates of /3(/x). 

We have, by (2.5), that 



/3(^) = ire-i^(-) 



ft / l-S 



eM^('''")&(q^)dr^ ds 






CM ^'^' 'b{qif)dr 



KQs) jo 



e 1^ 



^^(M>t) rt , 



A(q^) 

^A(m,s)' 

Qs). 



+ / e M 
s=0 -^0 



d / eM^^^'''^b(g^)dr- 



1 



/" eM^(^'^)6(q^)ds+ f 
Jo Jo 



* b(gn 



A(gr) ^^ 

Jo ^[Qs ) 

It is easy to see that 



X{q 



ds + e 



,-l7^(M>«) 



£.i*-.(,n*)^(^ 



|-R/3(M)llRd < 



|b||oo /" -hlU-s)j WHoo M 

^0 Jo ^0 ^0 



We also have 



(//) = _ / e-i^(^'^) ( / e-^^^''^'h{q>^)dr ] __VA(q^) • p^d. 



--Ai^J.,s) 



yo ^ J ^ kQs) ^ 

.e-i^(^'«) (p+- r eM^(^''-)6(q^)dr + - T e^^^^'^'UWr] ds 
\ fJ-Jo f^ Jo J 

(Ih) + (7/2) + (7/3) . 



Here 



(Ih 



t ---A{fi,s) 
^ ' ' ' e7^^^'''''h{q'^)dr]VX{qi^)-pds 



A2(q^) 







{Ih 



t -^A{^^,s) 



1 r e ^ 



^ Jo A2(qS 



eM^(^'")&(qnrf^ ) VA(q^^ 



eM^(^''')b(q^)dr ) ds 



(/^: 



1 / e f' ^ ^ 



/^Jo A2(g^) 
We conclude that 



^e^^(^'05(g/.)^^^ VA(g^) • ( T e'^^^^^^^Uw ^ ds 



\{IhMd < — -2 — \p\Rd\\b 

n n 9 

. ||VA||oo | /u^ 



loo ; e M ( / e m dr ) ds 



1 IIVAI 



|(//2)|m. <- ,2 

/i A, 



ll"lloo 



< 



IVA 







t / rs 



\J0 



{s-r)A, 



ii 



e 1^ dr \ ds 



00 III ||2 /^^ 
T2 ll"llooT2' ' 



E|(773)lL< f-^^llfolU) E 



1 IIVAI 



< 



< 



< 



< 



/^ 


Ag ' 


1 


|VA|U, 


/^ 


Ag ' 


1 


|VA|U| 


7t 


Ag ' 


1 


|VA|U, 


7t 


Ag ' 


l|VA||oo„, 
^ — 



/i A^ 

2 

&II00 1 E 



I V 

loo I 

t / fS 



\--.^m( Te-.^^^^'^^dr 



e-^^^^'''UWr ) ds 



g-i(A(M,s)-A(M,r)) 



dr I I / e f* 



V^o 

foiioo i'e fr (y"e-^(-^(^'^)"^(-''-))dryds] (f 



2-{A{f^,s)-A(p,r)) 



e 1^ 



dWr 



ds 







t / rs (s-r)xo ^ 2 



E 







g-^(A(M,s)-A(M,r)) 







t / rs 



e 1^ dr ] ds 
\Jo 

'^ I l-t / l-S (s-r)Ao 

e ** dr j ds I I / ( / e ^^ dr ) ds 

2 / , \ 2 / , 
t \ I lit 



dWr 



ds 



2(s-r)Ao 



An / V 2Ai 



Combining these estimates we see that E|(//)|^^ — )• as ^ J, 0. This imphes that 



E 



/3(/^) 



Ml 

A(qn 



ds 



as /Li 4, 0. 



2.3. Estimates of 7(^) - the reason why the classical Smoluchowski- 
Kramers approximation does not w^ork. 



We will show that E 



7(m) 



1 



Ata^""'- 



in general, does not tend to as 



'0 ^vys, 
/i 4, 0. Therefore the Smoluchowski-Kramers approximation does not work in the case 

of purely white noise perturbation. 

We have, by (2.5), that 



7(^) = i^e-^-^(-) 



e^^^^'^'UWr ] ds 



^ 



A(q: 



,M, / die-'.^^'^-^) 







A(qr 



-e '^ 



hAt^,t) 



Jo KKQs) Jo 



gM 



'^^^^''Uws . ,i 






-e '^ ^f"' ' _)_ 







gM 







A(q^ 



-dVrs+ 



dWr \d 







^7(^) + / TTZ^^dWs + (///) 



A(q: 



70 



/o A(q: 
It is easy to check that 



1 /-^ 2Xo(t-s) _ ^ 



^ ^0 



B\R^{fi%, <^ e M ds < -^ 



We have 



2A3 



Jo \Jo 

= {Ilh) + [Ilh) + {Ilh 



ei^(^'-)dW, 



X\q 



VXiq^)-p^ds 



where 



{Ilh 



t e-t^(M,s) 



A2(gS 



ei^('^''^)(iVl^ lVA(q^)-p(is, 



(Ilh) 



1 rt --A(u.,s) 
1 / e f* 



Ai io A^(q^ 



e^^(^'0^iy 1 VA(q^^ 



eM^(''''')6(qn^^ 1 (^s 



(^^^3) 



1 / e f* ^ ' 



We can estimate 



eM^^^-'^^dW, ) VA(q^^ 



e-.^'^f^'^UWr ] ds , 



B\iIII,)\l. < ' l^l-^ll^^l 



\2 
^0 



E 



--A{^l,s) 
e ^ 







e ** 



i{A{^,s)-A(M,r)) 



dl^r ds 



< 



< 



< 



|p|Rd||VA||c 

|p|Rd||VA||c 
-^0 



2 / --t _2Aji£ \ / ft / fS 

e 1^ ds 





jj_\ / fit 
2\J V2Ao 



vo 



g-l(A(M,s)-A(M,r))^ 



2^0(s-'-) 

e ^ dr 1 ds 



ds 



The term {III2) could be estimated in the same way as {Us)'- 

IVAI 



E|(///2)I4. < 



Aq J \^oJ V^-^o 



But in general one cannot estimate E|(///3)p up to a term which goes to as 
/i 4, 0. As an example, let A = ||A||oo and let us suppose that for 0<t<T<oowe 



have VA(q^) = ei. Here ei is the unit basis vector ei = (1,0, ...,0) in R.'^. Let WJ!^ be 
the fc-th (1 < A; < d) component of the Wiener process Wr- We have, for < t < T: 



E|(///3) 

1 



> 



yuA- 



E 



/iA2 



E 



\Jo 



e~M(^('^'^)-^(^''-))dW, 



g-l(A(M,.)-A(M,r))^^l , ^^ 



^-l(A(M,s)-A(^,r)) 1 , ^ . ^ 



+ E 



\Jo 



> 



1 



fc=2 
t / rs 



\J0 



^-l(A(^,.)^A(M,r))^^^fe 



1 



E 



^-i(A(;.,s)-A(M,r))^^^l , ^^ 



\^0 



^^-A^AM-AM)^^ ) ds 



(1 — e f* )ds 



;uA2 2A io 



M_ 

2A3 4A4 



^-i(A(M,s)-A(^,r))^ 1 , ^^ 



1-e M ) , 



10 



which does not tend to as /U J, 0. Since E|(///3)|^jj > (E|(///3)|jgd)^, we see that 
E|(I//3)|H^d does not go to as ^ 4, 0. Now we have 



E 



lip) 



HQs 



-dW. 



>-E\{IIh)\l,-E\R^{^l)\l,-E\{IIh)\l,-B\{III2)\i 



Therefore E 



7(^) 



constant as fi ^0. 



HQs 



-dW. 



is uniformly bounded from below by a positive 



We can check now that the process qf , < t < T, does not converge as fi I to 
the process Qt, Qo = Q- We have 



qr = g + 






ds + 



-dWs+ 



HQs)^ 



A(g^ 



-dW. 



/o -^(9j Jo HQs 



■dW.. . 



Suppose that we have, for any k,T > and any Pg = p G K fixed, that 

limP ( max Iqf — oJnd > k | = . 
We have, for some A> independent of /x and k, that 



Jo 
< AE max jq^ — Qg 



b(gn b(q. 



A(g^) A(qJ 



ds 



1 



VA(gs) A(qg 



dW. 



0<s<t 



<A 

yO<S<t ' 

< A[k + o(^,k)] , 



P max |q^ - qjL > k • E max |g^ - qjL + P max |q^ - qjL < k • k 



0<s<t 



0<s<t 



since E max |q^ — qr^|^^ < 00. Here the term o(/i, k) converges to as /i | for every 
fixed K > 0. Fix k > 0, let n l 0, we see that 



limE 

^lio 



[Qt - Qt) 



Kqt) b{q. 



vA(q^) A(gJ 
Since k > is arbitrary, we see that 



ds 



1 



A(q^) A(qJ 



dW. 



< Ak . 
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limE 

filO 



Jo 



7 5(9 






ds 



1 



1 



A(g^) A(g, 



dW. 



. 



On the other hand, let us suppose that VA(q^) = ei for < t < T < oo. Here ei 
is the unit basis vector ei = (1,0, ...,0) in W^. We have 



E 



°''''H«"'-i'ii'")+(^''"-rA(^)^"'' 



>1e 

- 3 



lip-) 



A(gS 



-dW. 



E|a(/i)|j 



E 



/3(/") 



A(gS 



-ds 



It follows from our estimates that this leads to a contradiction. 



3 Regularization via approximation of the Wiener process 

We could regularize the problem via approximation of the Wiener process. To this 
end we introduce the process 



1 



"^'-d,, '^•^ s 



s -t 



ds = \j^ W^^.pQds, 



where p(») is a smooth C°^ function whose support is contained in the interval [0, 1] 
such that 



p{s)ds = 1 



One can prove that (see [3] and the references there) 



limE max m^-WAl, 



. 



We have 



Wt 



Wt+5rP{'r)dr 



We can then introduce the following regularization of our problem: first we consider 
the system 



pq^/ = h{qf) - \{q^/)q^/ + wf , qf = q G 



l^ , qf =p(^W^. 



(3.1) 



Equivalently it is the first order system 



12 



• 11,5 u,5 

Qt =Pt 



pr = -M^r'') - ^^pr- + -w, 



■ Mi"? 



/^ 






(3.2) 



We can proceed with the estimates similar to in the previous sections. Since for 
fixed 6 > , 



\Wt\^d < - I m&K \p{r)\ ) ( max jVr^ljjd ) < oo , a. s. 

\0<r<l / \t<s<t+& 



(3.3) 



we could prove that all the terms 



E|a(/i)|Rd , E 



/3(/u) 



'h{q>i 



A(q^ 



-ds 



E 



7(^) 



t 1 

Kq'sf'^' 



goes to zero as /i J, 0. (To be precise, we should write a{fi,5), P{fi,5) and 7(;U,(5) to 
indicate the dependence on 5, but for brevity we neglect that.) In particular, with 6 > 
fixed, we can estimate the term {III3) up to a term which tends to as ^ J, 0. We have 



E|(//73)|m^ < -^^ / E 



1^ ^0 -^0 



^-L^AM-AM)^S 







ds 



1 IIVAI 
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Ai Ap ./n 52 



^0 ^0 



1 IIVAI 



/i Ap ./n 5^ 



E 



E 



^0 ^0 

1 ||VA||oo '■* 



\ ^{A{^--)-^(M,r)) ( f' Wr+5mP(.m)dm ] dr 











ds 



1 />s 

p{'m)Wr+5mdm / e 

JO 
/ -^ ( max |/5(m)| I E 



.i(A(M,s)-A(^,r)) 



dr 



2 / \ 2 

max iVt^/liR 
2 / \ 2 



ds 



^0 
<H TT^T^- ( max |p(m)| ) E( max IW; 

An 5^ V0<m<l / V0<Ks+<5 



e 1^ dr ] ds 



Therefore, for fixed 5 > 0, we have E| (///■? 



as /iiO. By (2.4), we get: 



H,<5 I 

Qt =1 + 



b{q 



fj.,S\ 



A(q: 



^,5 



-ds + 



1 



\{qn 



fi,5 



-dWi + 



+a{p) + Pin) 



' b{q^'') 



ds] + \ j{fi) 



fo KQs' J / 
Let the process q^ be governed by the equation 

,S _ b{ql) , 1 



A(g; 



-^dWl 



Qt 



+ 77^ ^t , q^ = g G M'^ 



Kq!) Kq! 



Then 



(3.4) 



(3.5) 
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'J 



A(qf) Jo \{q' 



Qt = <l+ T^ds+ / —rrdWt • (3.6) 



Let M{t,6,^) = E max \qs' — qflRd- By (3.4) and (3.6), using estimate (3.3), we 

0<s<t 

have 

M{t,5,ii) <Ki I M{s,6,fi)ds + K2it,6) [ M{s,6,fi)ds + 0^,(1) . 
Jo Jo 

Here o^(l) is a term which goes to as /i J, 0. The positive constant Ki is inde- 
pendent of /i, 6 and t. The positive constant K2 = K2{t, 6) may depend on t and 6, but 
is independent of /i. Now we use the Bennian-Gronwall inequaUty: 

M{t,6,fi)<o^{l)eMiKi + K2{t,6))t) . 
We conclude that for any 6, k,T > fixed and any Pq' = p fixed, 

hmP ( max Iq'^' — OiLd > «; | = . 
m40 \o<t<T' ' ^"^ J 

Now we can take 5 4, 0. Using Theorem 6.7.2 from [10] we get the following result. 

Theorem 3.1. We have, as 5 I 0, that 

liniE max Jg^-qjjjd =0 , 

where q^ is the solution of the problem 

Here the stochastic term is understood in the Stratonovich sense. 
In the general case 

^gf'^ = b{q^'') - X{q^'')q^'' + a{qf)wt , qf = q , q^;'' = p , (3.7) 

where the matrix a{») satisfy assumptions made in Section 1, we have, similarly, that 
for any 6,k,T > fixed and any pg' = p fixed, 

limP I max \q^' — q^Ld > k ) =0 
M^o Vo<t<T' ' ^"^ ' 

The process q^ is governed by the equation 
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7it='^ + ^wl4 = geR'. (3.8) 

A(g^) X(qt) 



And we conclude with 



Theorem 3.2. Under the assumptions mentioned above, 

lim E max \q^ — qJm,d = , 
5^0 te[o,T] 

where q^ is the solution of the problem 



^-li-ly"""*-^'''- <=■'' 



4 One dimensional case 

In the case of one space variable, Snioluchowski-Kraniers approximation leads to 
an one-dimensional diffusion process qt which is defined by the following stochastic 
differential equation written in the ltd form: 

<it='^-:^ + ^m,qo = ,eRK (4.1) 



X{qt) 2XHqt) X{qt) 



Put 



u^q) = / Hx) exp ( -2 / h{y)X{y)dy ) dx , 

v{q) = 2 r X{x) exp h J b{y)X{y)dy\ dx . 

Since A(x) > 0, u{q) and v{q) are strictly increasing functions. Following [4] we 
introduce an operator DyD^, where D^ means the differentiation with respect to the 

monotone function u(q): Duf(q) = lim —. — --—: the operator D^ is defined in 

?i->o u{x + h) — u{x) 

a similar way. One can check that D^D^ is the generator of the diffusion process qt 
defined by (4.1). 

Suppose now that the friction coefficient X{q) = Xs{q) depends on a parameter 
e > 0. We assume that, for each e G (0, 1], Xgi^q) has a bounded continuous derivative 
X'^{q), and < A < Ae(g) < A < oo. Let Ue{q) and Vir{q) be the functions defined by 
(4.2) when X{q) is replaced by Xe{q). 

Consider the stochastic process q^' '^ in M^ defined by the equation 

■■u,S,e 1/ u,5,e\ \e/ u,5,e\ u,S,e , jirS A'i'5,e •A'i<5,e i" a o\ 

m = Kit ) - A (it m +Wt ,10 =q , Qo =p ■ (4.3) 
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where W^ is, as before, a "smoothed" white noise converging to Wt as 6 ], 0. 

Theorem 4.1. Assume that the function Xs{q) converge weakly as e ^ on each 
finite interval [a, /3] C M^ to a function X{q) (maybe, discontinuous). Then processes 
q^' '^ converge weakly on each finite time interval to the diffusion process q^ governed by 
the generator D^Du (where u{q) and v{q) defined by (4-2) with A = X{q)) as, first /U | 0, 
then 5 10, and then e J, 0. 

Proof. According to Section 3, processes q^' '^ converge weakly as first /J, ], and 
then (5 4, to the process q^ which solves equation (4.1) with X{q) = X^{q). It follows 
from our assumptions that functions Us{q) and ^^(g) converge as e J, to functions u{q) 
and v{q) respectively for each g G M^. The functions u{q) and v{q) are continuous and 
strictly increasing. Therefore ([4]) a diffusion process q^ exists governed by D^D^. As 
shown in [8], convergence of Us{q) and V£{q) as e J, to u{q) and v{q) respectively implies 
weak convergence of processes qf to the process corresponding to DjjDu as e 4, 0. D 

Theorem 4.2. Let Xs{q) = A ( - I . Assume that one of the following conditions 
is satisfied: 

1. X{q) is a continuously differentiable positive 1-periodic function; 

2. X{q) is an ergodic stationary process (independent of the process Wt in (4-3)) 
with continuously differentiable trajectories and < A_ < A((7) < A+ < 00 for some 
constants A_, A-|-. 

Put A = / X{q)dq if condition 1 is satisfied, and X = EA(g) is condition 2 is 
Jo 
satisfied. 

Then the process q^' '^ defined by (5.3) converge weakly when first ^ ], and then 

e 4, to the process q^ defined by the equation 

Qt = jHQt) + jWt , qo = q ■ 



Proof of this theorem follows from Theorem 4.1 since each of conditions 1 and 2 
implies conditions of Theorem 4.1 and X{q) = A. D 

Assume now that Xs{q) is a bounded and separated from zero uniformly in e € (0, 1] 
positive function such that limAe(g') = Ai for q < 0, and limA£(g) = A2 for g' > 0. 

Assume that Xi;{q) is continuously differentiable for each e > 0. Let X{q) be the step 
function equal to Ai for g' < and to A2 for q > 0. Let functions u{q) and v{q) 
be defined by formula (4.2) with X{q) = X{q); u{q) and v{q) are continuous strictly 
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increasing functions. Denote by qt the diffusion process in M^ governed by the generator 

A = D^Du- The process qt behaves as —Wt on the negative part of axis q and as —Wt 

Ai A2 

on the positive part. Its behavior at q = is defined by the domain of definition S^ 

of the generator A: a continuous bounded function /(g), q € M-*^, twice continuously 

differentiable at g G {M^ \{q = 0}} belongs to I^a if and only if left and right derivatives 

at Q' = 0, /i(0) and /+(0) respectively, satisfy the equality — /i(0) = — /^(O) and 

Ai A2 

Af(q) is continuous. 

It is easy to see that functions Ue{q) and Ve{q) defined by (4.2) with X(q) = Xe{q) 

converge as e J, to u{q) and v{q) respectively for each g G M"^. This implies the following 

result. 



Theorem 4.3. Let the friction coefficient Xe{q) satisfies the conditions mentioned 
above. Then the stochastic process q^' '^ defined by (4-3) converges weakly to the diffusion 
process qt in E} governed by A = Dr^D^ as first fi iO, then 6^0, and then e 4, 0. 

This means, roughly speaking, that, if the friction coefficient is close to the step- 
function X{q), then process q^, for < fi « 1, can be approximated by the diffusion 
process %. 



5 Multidimensional case 

In this section we consider the problem of fast oscillating periodic environment in 
multidimensional case. We consider the system 

Here as in Section 3 the process W^ is the approximation of the Wiener process in 
W^. We make the same assumptions about the functions A(») and fo(») as in Section 2. In 
addition we assume that the functions A(») and 6(») are 1-periodic, i.e. X{x+ek) = X{x) 
and b{x + e^) = b{x) for a; € M'^ and e^ = (0, 0, ..., l(fc-th coordinate), ..., 0), 1 < k < d. 
Under this assumption our system (5.1) could be regarded as a system on the d-torus 
j^d _ ^d^^d^ Yix e > 0, we can proceed as in Section 3 to see that first as /U 4, then as 
(5 4, the process q^' '^ converges in probability to the process qf subjected to 

b('l 



^* = /n^\ + /»e\ °Wt ,ql = q(^ 
Af^) A(^ 
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The above equation, written in the form of Ito integral, will be 



gf 



6(| 



1 



VA(^ 



2.,3f| 



+ 



1 



A(^ 



Wt ,qf) = q€ 



l)d 



(5.2) 



e J \ £ / \ £ 

The generator corresponding to (5.2) is the second order differential operator 



L^uix) 



X 



1 



VA 



X 




1 1 



Vu{x) + ^ ^Au(x) 



u 



(5.3) 



A ( - ) 2e ^3 

.£/ \e. 

Our goal is to study the homogenization properties of (5.3) for general multidimen- 
sional case. Homogenization problems are considered by many authors, see, e.g., [7], 
[15], [14], [13], [12]. However, we provide here an elementary probabilistic way of doing 
this. Our method follows [7] and [6] (pp. 104-106). 

Let us first make a change of variable — = y and -^ = s. The process y^ = — qf 
corresponds to the generator 



A' 



1 



-A -XM.v +eM.v 

Ky) 



2A2(y) y 2XHy) ^ 



We regard y^ as a process on T . Then we have the bound 



^q/ef{yl) - I f{x)^f{x)dx 



<Ke 



Here K > Q and a > are independent of e for small e. The function / is bounded 
and measurable. The function fJ,^{x) is the density of the unique invariant measure of 
2/1 on T"* and / fjf{x)dx = 1. We have 

lim/i^(£c) = ;u(a;) , lim / f{x)fi'^{x)dx= / f{x)fi{x)dx 

elQ eiO Jjd Jjd 

for / G C(T^) and /i(£c) the unique invariant measure for the process with generator AP 
on T and / fi{x)dx = 1. Combining these estimates we have, that for any n, for any 
t > 6 > 0, there exist eo{n, 6) > such that for any < e < eo(^) <^)) we have 





^"f (f ) - I, f(.^)K^)d^ 


< - 


1 


any / G C(T'^), 


lim 


sup 
t>5 


E<,/ (^) - [ f{x)fx{x 


)dx 


= 



. 
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Finally we calculate the density fj,{x). Since 



^ A. - S4M • Vy = ^^(A^ - V(lnA(y)) • V^) 



2A2(y) y 2A3(y) 

we see that fj,{x) = CX{x) with C = [ X{x)dx I and we have the following result: 
Lemma 5.1. For any f € C{T'^), we have 



lim sup 



E„/ 



jd 



f{x)\{x)dx 



jd 



\{x)dx 



(5.4) 



Corollary. For any hounded continuous function f{x) on T^, q €z T'^ we have 

2 



E„ 



as e 10, for < t < oo. 



Ve 



t f{x)X{x)dx 



\{x)dx 







The proof of this corollary follows the same proof of the corollary after Lemma 1 



in [7]. 



Now let us consider auxiliary functions Nk{y), k = l,...,d, which are the periodic 
bounded solutions (i.e., on T'^) of the equations 



^ AyNkiy) - ^f^ • VyNkiy) = A^{Nk{y)) 



2X^{y) 



2\Hy) 



1 dX 

2X^{y)Wk 



:?/) , 2/ G T^ (5.5) 



The solvability of this equation comes from the fact that {AP)*X{y) = and 

1 dX 

{y)X{y)dy = 0. The boundedness of solution comes from our assump- 



jd 2A3(y) dyk 
tions on the function A(»). Now we apply Ito's formula: 
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.«. I f 



Ml 



ViVfe I ^ 1 • - 
\ e J e 



( 



V 



We 



2e A3 \e 



-t\ — \ -^^^\ — + 



W, 



ds+ 



1 

+ 2 



AiV,(^^ ^ 



e y £2^2 fl' 



r-(?) 



/ 



ds 



V 



Ale 



+ 



VT. 



A( ^ 



ds + 



2e 



dX 



J 



X3 1± 



ds 



(5.6) 



Let N{y) = (iVi(y), ...,Nd{y)). Using (5.5) we have 



/ 



(ft~q 



\ 






ds + / (DN) (^ 



b (q's\ , Ws 
aU /^^^Q| 



(is+ 



Here {DN){y) 



dNi 



y = iyi,...,yd)&T^. 



. '^Vj 7 l<i,j<d 

Therefore using the corollary after Lemma 1, we see that qf converges weakly to a 
process q^, Qq = q £ M."^ governed by the operator 



L 



d 



d^ A- d 

a.iA - — h / ^0 



Hj 



dyidyj 



«,j = l -^ z=l 



9t/i 



with coefficients 



(5.7) 



_ / VJV,(y)-ViV,(j/) 1 

"■ " ' V Ky) Ky) V % 



^«j 



fd 



dNj. , aiv. \ 1 



9yj 



-Jij 



Ky) 



dy 



jd 



Ky)dy 



/ bi{y)dy d / hk{y)^{y)dy 
Jfd sr-^ Jjd qyk 

\{y)dy k=i / \{y)dy 

Here 5ij = 1 if i = j, and 5ij = otherwise. 

We could simplify the expression for ojj: using (5.5) we get 



(5.8) 
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^ VAT.ta) ^ VNM ^ ^ (ON^^^^ ^ a«.j^j^ ^ ^.^,^ , ^^ / , ^ ^^^^^^ 



/ird V '^(y) Hy) \ 9yi dyj J \{y) J / \Jjd 

— iV, y -— - — -— ^ y + 
TTd \dyj V ^(Z/)/ A(y) a^j 

jd dyi X{y) Jjd X{y) 

~r 0- 



/ A(2/)dy '' / X{y)dy 



(5.9) 



So we have 



Theorem 5.1. As e \^Q, the process qf converges weakly to a process q^, qo = q & 
W^ governed by the operator (5.7) with coefficients given by (5.8) and (5.9). 

This Theorem imphes a homogenization result for the process q^' '*" defined by 
equation (5.f). 



6 Remarks and Generalizations 

6.1. Small mass - large friction asymptotics. 

Let the friction coefficient in (1.1) be A^(g) = e"^A(g), q G M", < e << 1. As 
it follows from Theorem 3.1, the Smoluchowski-Kramers approximation in this case has 
the form: 

Put ql = g w . Then q| satisfies the equation 



H/e- 



^* = A(i)-W(iy+A(i)'^*''°="' ^'-'^ 



where Wt is a Wiener process. 

Assume that the vector field b(q), q £ M", has a finite number of compact attractors 
Ki,...,Ki. Let, for brevity, each Ki be an asymptotically stable equilibrium, and each 
point of M", besides a separatrix set £ C M", is attracted to one of these equilibriums. 
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The separatrix set £ is assumed to have dimension less than n. Then, if Qq = q f, ql 
first comes to a small neighborhood of a stable equilibrium Ki, i = i{q), with the proba- 
bility close to 1 as e 4, and spends in this neighborhood a long time. Then, because of 
the large deviations, the trajectory will switch to the neighborhood of another attractor, 
then to another, and so on. We see from (6.1), that the long-time behavior of the system 
with a large friction is similar to the behavior of a system with a small noise. Applying 
the results of [9], Chapters 4, 6, we see that, for < e << 1, the sequence of transi- 
tions between the attractors, the main term of transition time logarithmic asymptotics, 
and the most probable transition paths are not random for generic systems. These 
characteristics of the long-time behavior are defined by a function V{x,y): 

V{x,y) = mfi^l \X{.^,)ip,-b{if,)\^ds:ipo = x , ifT = y,T>o\ , x,y eW ; 

and by the extremals of this variational problem. 

Let now the dynamical system q^ = b{q^) has a first integral. Let, say, n = 2 and 
6(g) = \7H{q) for some smooth generic function H{q), q € M?, such that lim H{q) = 

\q\->-co 

00; in this case H{q^) = H{qQ). 

Assume again that the friction is strong: A^(q) = e~^X{q). Make a time change 
Qt =Qt/e'^- Then 

Identify points of each connected component of every level set of H{q). The set 
obtained after such an identification if homeomorphic in the natural topology to a graph 
r. Let y : M^ — 7- r be the identification mapping. Then the long-time evolution of the 
system can be characterized by the stochastic process yf = Y{ql) on F. The process 
yf, in general, is not Markovian. But yf converges weakly in the space of continuous 
functions (/9 : [0,T] — )• F as e J, to a diffusion process on the graph F ([9], Ch.8). This 
limiting process is defined by a family of second-order differential operators, one on each 
edge of F, and by the gluing conditions at the vertices. Following [9], one can evaluate 
these operators and the gluing conditions. 

6.2. Fast oscillating random friction in multidimensional case. 

Let us consider the case of fast oscillating in the space variable, random friction, in 
dimension d > 2. Let (0,, J-", P) be a probability space. Let A(a;, oj), u £ Q he a random 
field in M with the following properties: 

(i) For any fixed a; G fi and x G W^, the function 00 > A > X{x, 00) > Xq > 0. 

(ii) For every x (^ M.'^ the random variable X{x,uj) is independent of the Wiener 
process Wt- 
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(iii) The random field X{x, lo) has the form A(a;, cj) = X{T{x)uj) where T{x) : il — t- $7 
is a d-dimensional dynamical system which preserves the measure P and is ergodic with 
respect to P. 

Let us now consider an analogue of (5.1): 



■■ u,5,e 



Qi 



IJ,,S,e 



u qf'^ + W\ , qf^' = q G M^ q^/^' = p G M'^ 



3.2) 



For each fixed a; G fi, as we proved in Section 3, we have that qf' '"^(w) converges 
weakly to a process qf (w) as first ^ | and then 5 | 0. The process qf is subject to 



gf 



1 



VA( ^,u: 



''X^('§,u 



+ 



Xi^.u 



Wt , qo = q G 



(6.3) 



We conjecture that as e J, 0, qf converges weakly to a process q^, qo = q G 
= 1 d ^ Q2 
subject to the operator -^ = - ^ ^ij~ — KZ~ with effective diffusivity 



«J=i 



dxidxj 



a-ij — E 






[X,UJ 



X{x,c 



-dx 



+ <5, 



fd X{x^oj 



-dx 



*j 



(6.4) 



X{x,ijj)dx I X{x,ijj)dx 

Ijd Jjd 

Here the functions Nk{x^uj) {1 < k < d) shall satisfy certain auxiliary problem. 
(We actually have a formulation of this problem but we are not sure about its validity: 
we let Nk{x,uj) be the solution of the equation 



^[{VxNk{x, w) - Sk) ■ Vx'fiix, tj)] = , 

for all (p{x,co) smooth and compactly supported in a; G M and measurable with respect 
to a; G f^. The existence of solutions to this problem is guaranteed by the Lax-Milgram 
lemma.) 

However, we are not aware of the validity of this conjecture nor a proof of it. We 
are also not sure about the correct reference of such a problem. (We thank E.Kosygina 
for pointing out to us two relevant papers [1] and [11].) 

6.3. Motion of charged particles in a magnetic field 

One can expect that using the regularization by smoothed white noise one can get 
the Smoluchowski-Kramers approximation for the equation 



Mi = Hqi) - MQt)(it + o-W't , Qo = g , 90 = p ; q,p e 



(6.5) 
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Here a > and A{q) are a matrix-valued functions having strictly positive eigen- 
values for each q £ M". In particular, if A{q) = X{q)A, where A(g) > and A is a 
constant positive definite symmetric matrix, the problem can be reduced to the case 
considered in Section 3 by an appropriate linear change of variables. 

If A has a negative eigenvalue, the Smoluchowski-Kramers approximation is not 
applicable. The case of A with pure imaginary eigenvalues is of interest since such 
equations describe the motion of charged particles in a magnetic field. If A = const, n = 2 

and ^ = , the problem was considered in [3] . In this case the Smoluchowski- 

Kramers approximation holds after a regularization. If 6(q) = —\7F{q), q € M?, and 

. (^ -A 

A = 

to the equation 



A = , one can show that the regularization by the smoothed white noise leads 

\1 / 



If the noise in (6.5) is small (0 < o" << 1) the motion described by (6.6) has a fast 
and a slow components. Applying the results of [9], [2], one can describe the limiting (as 
a I 0) slow component of q1 = qt/„2 as a diffusion process on the graph corresponding 
to the potential F{q) (the graph is homeomorphic to the set of connected components 
of the level sets of F{q) provided with the natural topology). 

Acknowledgements: This work is supported in part by NSF Grants DMS-0803287 
and DMS-0854982. 



References 

[1] G.Allaire, R.Orive, Homogenization of periodic non self-adjoint problems with 
large drift and potential, Esaim-control Optimisation and Calculus of Variations , vol. 
13, no. 4, pp. 735-749, 2007. 

[2] M.Brin, M.Preidlin, On stochastic behavior of perturbed Hamiltonian systems, 
Ergodic Theory and Dynamical Systems, 20 (2000), No.l. 55-76. 

[3] S.Cerrai, M.Freidlin, Small mass asymptotics of a charged particle in a mag- 
netic field and long-time infiuence of small perturbations. Journal of Statistical Physics, 
submitted. 

[4] W. Feller, Generalized second order differential operators and their lateral con- 
ditions, III J. Math., 1 (1957) pp.459-504. 

[5] M.Preidlin, Some Remarks on the Smoluchowski-Kramers Approximation, Jour- 
nal of Statistical Physics, 117, No.314, pp. 617-634, 2004. 

24 



[6] M.Freidlin, Functional Integration and Partial Differential Equations, Princeton 
University Press, 1985. 

[7] M.Freidlin, Dirichlet's problem for an equation with periodic coefficients de- 
pending on a small parameter. Theory of Probability and Its Applications, 9, pp. 133-139, 
1964. 

[8] M.Freidlin, A.Wentzell, Necessary and Sufficient Conditions for Weak Conver- 
gence of One-Dimensional Markov Processes. The Dynkin Festschrift: Markov Processes 
and their Applications , Birkhauser, (1994) pp. 95-109. 

[9] M.Freidlin, A.Wentzell, Random perturbations of dynamical systems. Springer, 
1998. 

[10] N.Ikeda, S.Watanabe, Stochastic Differential Equations and Diffusion Pro- 
cesses, second edition, North-Holland Publishing Company, 1989. 

[11] M.Klepsyna, A.Piatnitski, Averaging of non-self adjoint parabolic equations 
with random evolution, Institut National de Recherche en Informatique et en Automa- 
tique, June 2000. 

[12] N.Krylov, G-convergence of elliptic operators in non-divergence form, Math 
Zametki, 37, 4, pp. 522-527, April, 1985. 

[13] G.Papanicolaou, S.R.S.Varadhan, Boundary value problems with rapidly os- 
cillating random coefficients, Seria Coll. Math. Soc. Jdnos Bolyai, North Holland, 
Amsterdam, 1979. 

[14] V.V.Zhikov, S.M.Kozlov, O.A.Oleinik, Homogenization of differential operators 
and integral functionals. Springer, 1994. 

[15] V.V.Zhikov, S.M.Kozlov, O.A.Oleinik, and Kha T'en Ngoan, Averaging and 
G-convergence of differential operators, Russian Math Surveys, 34, 5, pp. 65-133, 1979. 



25 



